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A recent hypothesis claims that the glass transition itself, though it is a very 
pronounced relaxation peak, is no separate relaxation process at all, but is 
just the breakdown of the shear modulus due to the weak elastic dipole inter- 
action between all the quasi-independent relaxation centers of the glass. Two 
derivations are considered, one of them in terms of a breakdown of the shear 
modulus and the second in terms of a divergence of the shear compliance. 
Mechanical relaxation data from the literature for vitreous silica, glycerol, 
polymethylmethacrylate and polystyrene are found to be consistent with the 
first hypothesis. 

PACS numbers: 64-70. Pf 

In 1972, Sigi Hunklingers famous experiment 1 demonstrated the two- 
level nature of the universal low-temperature excitations in glasses. 2 The 
experiment gave strong support to the energy landscape concept, 3 in which 
the two-level states are explained in terms of tunneling between adjacent 
minima of the energy landscape, with a low barrier between them. 4 

If the energy landscape concept is indeed the correct explanation, these 
two-level states are expected to be a small fraction of a vast number of 
local relaxation centers in the glass. In this view, the relaxation spectrum 
of the energy landscape spans a wide range, from the tunneling states at 
the low-barrier side to the Johari-Goldstein peak 5 at the high-barrier side. 
Depending on frequency, this Johari-Goldstein peak is either observed in the 
glass close to the glass transition temperature T g or even above T g in the 
supercooled liquid. 

While the energy landscape concept in itself seems to be reasonably 
well established, the extent to which it can be represented by an ensemble 
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of independent double-well potentials is a widely debated question. With 
respect to the tunneling states, this question begins to be answered by ex- 
periments in the mK range. 6 Again, Sigi Hunklinger played a central role 
in the performance and interpretation of these experiments, approximately 
thirty years after establishing his scientific reputation with the demonstra- 
tion of two- level states. 1 According to these new experiments, the concept of 
isolated tunneling states seems to break down in the mK range. A plausible 
explanation is the coupling of the tunneling states by the interaction of their 
elastic dipole moments, 7 which leads to the formation of coupled pairs of 
tunneling states. 

The present paper deals with the question whether the same elastic 
dipole interaction provides the key to the unsolved riddle of the glass transi- 
tion. Here, we assume that this is indeed the case. In particular, we assume 
that the shear relaxation is due to an ensemble of double-well potentials 
which are weakly coupled by the interaction of their elastic dipole moments. 
The ensemble comprises the tunneling states at the low-barrier end as well as 
the Johari-Goldstein relaxation at the high-barrier end, possibly with rather 
different atomic jump vectors, but all of them interacting with each other. 
The reader should be aware that this is only one of many explanations in 
the literature. 8,9 

One can treat the elastic-dipole interaction in two different mean-field 
approaches (i) in terms of the shear modulus (ii) in terms of the shear com- 
pliance. Here, equations for both approaches are derived and compared to 
dynamical mechanical shear data for four glassformers at their respective 
glass temperatures T g . 

In the shear modulus treatment, 10 one considers a small initial shear 
deformation e switched on at the time t = 0. The initial shear stress a = Ge, 
where G is the infinite frequency shear modulus. The shear stress decay with 
increasing time is described 11 in terms of the rheological function H{t) 



In the energy landscape concept, relaxation occurs via thermally ac- 
tivated jumps over the barriers V between different energy minima. The 
corresponding relaxation time ry obeys the Arrhenius relation 



with To ~ 10 13 s. Therefore we replace the rheological function H(r) by the 
barrier density function f(V) 




(1) 



t v = r exp(V/k B T) 



(2) 



H( TQ e v / kBT 



) = H(t v ) = Gk B Tf{V). 



(3) 
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Inserting this definition into equation (1) for the time dependence of 
the shear modulus gives 

G(t) = G (l - J™ f(V)(l - e- t '^)dV S j « G (l - £ f(V)dvj , (4) 

where Vt is the barrier with ry = t. 

In order to take the influence of the elastic dipole interaction into ac- 
count, one assumes a true barrier density foiV) of the energy landscape. 
At the barrier height Vf, the jump between the two minima will tend to 
occur at the time t, when the shear stress has decreased to G(t)e. Thus 
the remaining shear stress energy is only a fraction G(t) 2 /G 2 of the initial 
one. If we now consider the microscopic situation at the local relaxation 
center as unchanged, with the same initial shear distortion, then the free en- 
ergy change by the jumps of the relaxation center remains unchanged. This 
means that the relaxation center releases the same amount of stress energy. 
But it reduces a stress energy weakened by G(t) 2 /G 2 , so its cffectivity in 
bringing the shear modulus down to zero is increased by the reverse of this 
factor. 

Therefore the mean-field approach for the shear modulus reads 



fo(V) = f(V) 



- / f(v)di 
Jo 



(5) 



In a physical picture, the enhancement of fo(V) is due to induced jumps 
of lower-barrier relaxation centers in the neighborhood of the given relax- 
ation center with barrier V, which occur quasi-instantaneously after its jump. 

The equation (5) has the back-transformation 10 

f( V ) = 77^ —y (6) 



1-3.^ fo(v)dv 



In the shear modulus approach, the breakdown of the shear rigidity 
occurs when the integral of fo(V) over V reaches 1/3. The corresponding 
barrier is called Maxwell barrier, because it determines the Maxwell time tm 
(the shear stress relaxation time) through the Arrhenius relation, eq. (2). It 
is given by the 1/3-rule derived in ref. 10 

V M 1 

fo(V)dV = -. (7) 



The breakdown occurs in a rather dramatic way, because the relaxing 
entities at the critical Maxwell barrier value receive a strong enhancement, 
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to such an extent that one is tempted to assume a separate a-process which 
has nothing to do with the secondary glass relaxations. In fact, this more 
or less unconscious assumption underlies most of the present attempts to 
understand the glass transition. 8 The above treatment shows such an as- 
sumption to be unnecessary; what one sees at the glass transition are simple 
Arrhenius relaxations of no particularly large number density fo(V), blown 
up to impressive size in f(V) by the small denominator of eq. (6). 

The shear modulus approach can be tested by plotting the product 
G"G' 2 of experimental dynamical shear data as a function of frequency or 
temperature. In fact, the true barrier density /o(V) can be calculated ap- 
proximately from the equation 

V M f (k B Tln— ) « -^fln^, (8) 

where G is again the infinite frequency shear modulus. 

One can build the analogue of the shear modulus approach, starting 
from the shear compliance. In the shear compliance treatment, one considers 
a small initial shear stress a switched on at the time t = 0. The initial shear 
distortion e = Ja, where J = 1/G is the infinite frequency shear compliance. 
The shear distortion increase with increasing time is described 11 in terms of 
the rheological function L(t) 

J(t) = n+ L (r)(l ~ e-'/^lnr + -, (9) 

where r\ is the viscosity. In the following, we omit this viscosity term, because 
we are only interested in the breakdown of the shear rigidity. 

In the energy landscape concept, we replace the rheological function 
L(t) by another barrier density function l(V) defined by 

L(r e v ^ T )=L(ry) = k -^l. (10) 

Inserting this definition into equation (9) for the time dependence of 
the shear compliance gives 



roo 1 1 rVt 

1+ l(V){l-e- t / Tv )dV -77 1+/ l(V)dV 
Jo J G Jo 



(11) 



where again Vt is the barrier with ry = t. 

Again, the influence of the elastic dipole interaction is taken into account 
assuming a true barrier density Iq(V) of the energy landscape. But now, the 
situation at the barrier height Vt for the time t is different. The shear 
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stress is still the initial one, but the shear distortion has increased by J(t)G. 
Therefore the free energy change by the jumps of the relaxation center is 
increased by J(t) 2 G 2 . 

Thus the mean-field approach for the shear compliance reads 



l (V) EE 

Its back-transformation is 
l(V) = 



l(V) 



l + J V l(v)dv 
k(V) 



2 ' 



(12) 



1 - Jo h{v)dv 



2 ' 



(13) 



We see that the results for the shear compliance approach (ii) differ 
markedly from the shear modulus approach (i) proposed earlier. 10 The 
Maxwell barrier is now given by the divergence of the compliance 



rV M 
Jo 



l (V)dV 



(14) 



Again, the relaxing entities at the critical Maxwell barrier value receive 
a strong enhancement. The shear compliance approach can be tested by 
plotting the ratio J"/ J' 2 of experimental dynamical shear data as a function 
of frequency or temperature. The true barrier density ^o(^) can be calculated 
approximately from the equation 



V M lo{k B T\n—) « - J" J /J' 2 In—. 



15 
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Since the influence of the viscosity was neglected, this equation is only ex- 
pected to hold down to the frequency with ujtm = 10, at least for non- 
polymeric glass formers. This is a drawback of the compliance approach. If 
one tries to include the viscosity, the equations get complicated and rather 
difficult to handle. 

In any case, one is now able to compare dynamical mechanical shear 
data from experiment to the two theoretical approaches. Fig. 1 begins 
with the shear modulus approach, showing G"'G" 2 -values from measurements 
at the glass transition. The comparison is done for the four glass formers 
Si02, glycerol, polystyrene and polymethylmethacrylate. Also included is 
the simplest theoretical case, with fo(V) = constant = 1/3Vm, the so-called 
generic case. The figure is scaled in such a way that this generic case is 1. 
The scaling requires the knowledge of the Maxwell time tm- For silica and 
glycerol, the classical Maxwell relation tm = rj/G was used. For the two 
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polymers, t\j was taken from the condition uj max TM = 1, where uj max was 
the maximum of the peak in G" . In order to cover five decades in such a 
mechanical measurement, one usually needs to do measurements at several 
temperatures around T g and use the time-temperature scaling to obtain the 
full curve. 

The scaling G"G'/G S also requires the knowledge of the infinite fre- 
quency shear modulus G. Therefore one needs not only a measurement of 
G' and G", but also a Brillouin measurement of the transverse sound wave 
velocity vt at T g , to determine G via G = pv%, where p is the density at T g . 
Table I summarizes the T g and G values of the four glass formers, together 
with the relevant references. The glass temperature for glycerol is somewhat 
higher than the usual value, because the master curves for G' and G" were 
given 17 for this temperature. 



Table I: Glass transition data. 



substance 


Si0 2 


glycerol 


polystyrene 


polymethylmethacrylate 


Tg (K) 


1449 


192.5 


363 




383 


G{T g ) (GPa) 
V (T g ) (GPa s) 


33.8 


4.3 


1.5 




1.86 


2512 


10 








TM (s) 


74.3 


2.3 


100 




100 


ref. G 


12 


13 


14 


15 




ref. G', G" 


16 


17 


18 


19 




m 


20 


53 


138 




145 



In Fig. 1, all four experimental curves show a cutoff at ujtm = 1. If one 
looks closely, one observes that this experimental cutoff is slightly broader 
than in the theoretical calculation for the generic case, where a sharp cutoff 
of /o(V) at Vm was assumed. The two polymers lie decidedly lower than the 
generic case, the two other cases lie higher. 

The corresponding comparison to the compliance approach in Fig. 2 
suffers from the fact that for silica and glycerol, it is only meaningful above 
ut~m = 10, because of the neglected viscosity. Again, silica shows the 
strongest rise towards the value ujtm = 1, PMMA the weakest. 

But even for the two polymers, which have such a high viscosity that it 
can be indeed neglected, the compliance approach does not work well; there 
remains a pronounced peak in the near neighbourhood of ujtm = 1, much 
more pronounced than in the shear modulus approach of Fig. 1. Therefore 
we conclude that the shear modulus approach is the better description. 

This is not unexpected: The shear modulus approach corresponds to 
a generalized Maxwell model, with an infinite number of parallel Maxwell 
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Fig. 1. Comparison of data measured at the glass transition of S1O2, glyc- 
erol, polystyrene (PS) and polymethylmethacrylate (PMMA) to the generic 
case fo(V) = 1/3 Vm in the shear modulus approach. 

elements, while the shear compliance approach corresponds to an infinite 
number of Voigt elements in series. 11 The description of an ensemble of 
localized relaxation centers in a threedimensional viscoelastic continuum in 
terms of parallel elements is probably more suitable. 

Further support for this conclusion is supplied by the obvious connec- 
tion between the values of the fragility m in Table I (taken from ref. 20 ) 
and the values of Vm/o{Vm) extrapolated from Fig. 1. The fragility 
m = d log i]/d(Tg/T) is defined in terms of the steep rise of the viscosity 
r] towards the glass temperature T g with decreasing temperature in the su- 
percooled liquid. Silica has the lowest fragility of all known glass formers, 
while the two polymers belong to the very high fragility end of the scale. 
A high fragility means a strong decrease of the Maxwell barrier Vm with 
increasing temperature. Fig. 1 indicates that a high m seems to be related 
to a low Vm/o^Vm), and viceversa. 

From the 1/3-rule, eq. (7), such a behaviour is expected: If fo(V) tends 
to increase with increasing temperature, a low value of /o(Vm) means a fast 
decrease of Vm with increasing temperature. Thus one gets at least a qualita- 
tive understanding of the old fragility riddle: 8 With increasing temperature, 
more and more minima of the energy landscape get populated. This leads to 
an increase of fo(V), which in turn leads to a decrease of Vm- This decrease 
is more dramatic for a low value of Vm/o(Vm) than for a high one. 
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Fig. 2. Comparison of data measured at the glass transition of Si02, glyc- 
erol, polystyrene (PS) and polymethylmethacrylate (PMMA) in the shear 
compliance approach. 

Of course, Fig. 1 does not prove the shear modulus description of the 
glass transition 10 beyond any possible doubt. But the plot of Fig. 1 is a new, 
quantitative and probably meaningful way to study the relation between the 
primary relaxation (the shear stress relaxation to the value zero) and the 
secondary relaxations at higher frequencies. 
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